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ABSTRACT 


'1C  Is  shown  Chat  under  natural  economic  assumptions 
the  complementary  pivot  algorithm  will  find  an  op¬ 
timal  steady  state  for  a  von  Neumann  growth  model 
with  a  quadratic  utility  function.  This  generalizes 
the  restilt  of  DanCzig  and  Mann  for  the  case  of 
linear  utility.  . 


A  COMPLEMENTARITY  ALGORITHM  FOR  OPTIMAL  STATIONARY 
PROGRAMS  IN  GROWTH  MODELS  WITH  QUADRATIC  UTILITY 


by 

’  Chla-Shln  Chung  and  David  Gale 

1.  INTRODUCTION 

Some  years  ago  Sutherland  [5]  and  Koopmans  [3]  Independently 
proved  the  existence  of  optimal  stationary  programs  for  growth  models 
In  which  future  utility  Is  discounted  at  a  constant  rate.  They  also 
showed  that  under  suitable  assumptions  such  programs  could  be  characterized 
as  ones  which  satisfy  a  simple  one  period  "competitive  condition."  Sub- 
4  sequently,  Hansen  and  Koopmans  [2]  gave  another  proof  of  this  charac¬ 

terization  and  showed  how  techniques  for  calculating  approximate  fixed 
*  points  of  continuous  mappings  could  be  used  to  approximate  these 

stationary  solutions.  An  alternative  approach  was  then  proposed  by 
Dantzlg  and  Mann  [1]  who  showed  that  for  the  special  case  of  a  linear 
utility  function  one  could  obtain  the  optimal  stationary  program  by 
applying  Lemke's  complementary  pivot  algorithm.  They  then  proposed 
approximating  an  arbitrary  concave  utility  by  piecewise  linear  utilities 
to  handle  the  general  case. 

The  main  purpose  of  this  note  is  to  show  that  the  result  of  [1] 
remains  valid  when  the  linear  utility  function  Is  generalized  to  any 
concave  quadratic  function.  This  represents  an  improvement  from  the 
economic  point  of  view  since,  as  has  often  been  observed,  linear  functions 
are  unsatisfactory  for  representing  utilities  since  they  Imply,  at  least 
for  Leontlef  technologies,  that  only  one  good  will  be  consumed.  Quadratic 
functions,  on  the  other  hand,  may  be  thought  of  as  a  good  first  approxi¬ 
mation  to  a  general  concave  utility  function. 
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A  secondary  objective  Is  to  reformulate  the  "key  hypothesis"  of  [  1  ] 
in  a  form  in  which  Its  economic  meaning  becomes  more  evident.  Finally, 
our  proof  has  the  merit  of  brevity.  We  assume,  however,  familiarity 
with  the  complementary  pivot  algorltim  whose  essential  properties  will 
be  glvctn  In  Section  3. 
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2.  THE  MODEL  AND  THE  COMPETITIVE  CONDITIONS 

We  are  given  m  by  n  matrices  A  and  B  ,  6  >  0  ,  an  m-vector  b  , 
a  concave  differentiable  (utility)  function  u  defined  on  and  a  dis¬ 

count  factor  6  on  (0,1)  . 

A  pTKgrcan  is  a  sequence  (x^)  of  non-negative  n-vectors  such  that 

(2.1)  ^t+1  ~  ®*t  ^  ^  . 

The  value  v  of  the  program  (x^)  is  given  by 

00 

v((x  ))  ■  \  6*^u(x  )  . 

t-0 

%  A  program  is  stationary  if  x^  -  x  for  all  t  and  is  optimal 

stationary  if 

v((x))  >_  v((x^)) 

for  all  programs  x^  with  x^  >■  x  . 

Definition; 

The  stationary  program  (x)  is  competitive  if  there  is  a  non-negative 
(price)  m-vector  p  such  that 

(2.2)  (A  -  B)x  <  b 

and 


(2.2)’ 


p(A  -  B)x  •  pb 
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t 


(2.3)  p(A  -  6B)  >  7a»(x) 
and 

(2.3) '  p(A  -  (5B)x  *  x7u(x)  . 

The  term  coo^etltlve  reflects  the  fact  that  the  conditions  above  are 
satisfied  precisely  when  x  maximizes  "profit",  defined  as  (i)(x)  + 

5pBx  -  pAx  ,  among  all  x  satisfying  (2.2). 

The  following  basic  result  Is  also  proved  In  [5]  and  [2].  For  the 
sake  of  completeness  we  reprove  It  here. 

Theorem  1; 

A  competitive  program  Is  optimal. 

■Let  (x^)  be  any  program  with  x^  «  x  where  (x)  is  competitive. 
From  (2.3)  we  have 


x^V(o(x)  <_  p(A  -  6B)x^ 

and  multiplying  by  5^  and  summing  from  0  to  T  gives 


T  T 

\  6  X  7u(x)  ^  p  \  6^(A  -  6B)x 

t-0  t-0 

-  P^Ax^  +  \  6^(Ax^  -  Bx^_j^)  -  S^'^’^Bx^j 

(from  (2.2)  and  B  >  0) 

T 

■  pAx  +  p(A  -  B)x  1  6^  (from  (2.2)'  and  x  »  x) 

t-1  ° 

T 

-  p(A  -  6B)x  I  5^  + 

t-0 


(i.-) 


T+1 

x7u(x)  +  6^  pBx  (from  (2.3)') 
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so 

T 

I  5*^(x.  -  x)7(o(x)  <.  6'^'^Sbx 
t-0 

T+1 

aod  since  6  pBx  approaches  zero  we  have 

m 

I  6^(x  -  x)7uj(x)  <  0  . 
t-0 

By  concavity  u)(x^)  -  a>(x)  <  7u)(x)(x^  -  x)  so 

00  09 

I  6‘^a»(x  )  _<  I  5^(ii(x)  or  v(x)  >  v(x  ) 
t-0  ^  t-0 

so  (x)  Is  optimal  ■ 

The  converse  of  Theorem  1  does  not  hold  as  shown  by  the  following 
simple  example: 

A-l,B-3,  u(x)  -  x>b»0  , 

The  trivial  program  x  -  0  Is  clearly  optimal  since  It  Is  the  only 
program  starting  from  x^  -  0  ,  but  It  is  not  competitive  since 
A  -  6B  -  -y  ,  so  (2.3)  has  no  solution. 

To  get  around  this  type  of  problem,  Peleg  and  Ryder  [4]  have  defined  a 
technology  to  be  i-produative  If  It  contains  an  Input-output  pair  (x,y) 
such  that  X  <  6y  .  This  Is  a  rather  natural  condition  which  asserts 
that  the  model  Is  capable  of  expanding  at  a  rate  greater  than  the  reciprocal 
of  the  discount  race.  For  our  purposes  Che  weaker  assumption  x  <  dy  will 
suffice.  Intuitively,  If  this  condition  were  not  satisfied  It  would  mean 
that  Che  growth  potential  of  the  technology  was  not  sufficient  to  compensate 
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for  constimers'  preference  for  present  over  future  consumption.  Thus, 
capital  would  be  consumed  faster  than  It  could  grow  and  the  only  possible 
steady  state  would  be  at  the  subsistence  level.  (For  b  >  0  the  only 
solution  would  be  the  trivial  steady  state  x  •  0.) 

In  our  model  the  (S-productlvlty  condition  becomes 

(2.4)  Ax  -  b  <  6Bx  for  some  x  ^  0  . 

A  second  assumption  we  will  make  Is  that  the  model  Is  not  capable  of 
achieving  steady  states  with  arbitrarily  high  utility.  This  Is  certainly 
reasonable  for  any  economy  In  which  production  depends  on  exogenous 
bounded  resources.  Formally  we  have 

(2.5)  There  exists  M  >  0  such  that  a)(x)  for  all  x  satisfying  (2.2). 

Condition  (2.4)  and  (2.5)  will  be  used  here  Instead  of  the  "key 
hypothesis"  used  In  [1]. 
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3.  THE  CASE  OF  QUADRATIC  UTILITY 

We  now  consider  the  special  case  where  cd  Is  concave  quadratic, 
that  is, 

u(x)  *  cx  -  y  xDx 

where  D  Is  symmetric  and  positive  seml-deflnlte.  In  this  case 

7u(x)  >  c  -  Dx  .  Substituting  this  In  (2.3)  and  (2.3)'  and  Introducing 

slack  variables  u  *  b  -  (A  -  B)x  ,  v  ■  p(A  -  6B)  +  Dx  -  c  ,  we  get 

the  competitive  condition  In  the  form  of  the  following  Linear  Complementarity 

Problem: 

Find  non-negative  (m  +  n)-vectors  z  »  (x,p)  and  w  “  (u,v)  such 

that 

(3.1)  w-Mz+q  ,wz*0 

where  ^  ).  q-  {~l)  • 

Our  main  result  Is 

Theorem: 

If  (2.4)  and  (2.5)  are  satisfied,  then  (3.1)  has  a  solution. 

■Recall  that  the  complementary  pivot  algorithm  Introduces  the  auxiliary 
vector  e  all  of  whose  coordinates  are  one,  and  either  finds  a  solution  of 
(3.1)  or  else  non-negative  solutions  (w,z,6)  and  (w,z,6)  of 


(3.2) 

V 

<i> 

1 

1 

> 

■  q 

,  wz  -  0  , 

0 

>  0 

(3.3) 

w  -  Mz  -  0e 

-  0 

,  wz  ■  0  , 

z 

1*  0 

4 


(3.4) 


W2  ■  WZ  «  0  . 
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We  must  show  that  under  our  assumption  this  second  alternative 
cannot  occur.  The  key  fact  we  need  is  that  the  matrix  M  is  "co¬ 
positive",  that  is. 

Lemma: 

If  z  >  0  ,  then  zMz  >  0  . 

s  at 

■  M  +  j  and  -  6)B  is  non-negative. 

Hence  if  z  «  (x,p) 

2(zMz)  ■  z(M  +  M^)z  =  2(xDx  +  p(l  -  6)Bx)  >  0 
since  D  is  positive  seml-def inlteB 

Now  suppose  (3.2)  -  (3.4)  hold.  Multiplying  (3.3)  by  z  gives 

zMz  +  0(ez)  *  0 

so  from  the  lemma,  we  get  9=0  and  zMz  *  0  since  z  0  •  Letting 
z  =  (x,p)  this  implies  xDx  *  0  which  means 

(3.5)  xD  =  0 

by  the  property  of  positive  semi-definite  matrices.  Since  9=0, 
(3.3)  becomes 

(3.6)  w  =  Mz 

so  multiplying  by  z  and  using  (3.4)  we  get 

(3.7)  zM^z  -  0  . 
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Multiplying  (3.2)  by  z  gives 

e 

(3.8)  -zMz  -  9ez  ■  qz 
and  adding  (3.7)  and  (3.8)  gives 

(3.9)  -z(M  +  M^)z  -  6(ez)  •  qz  . 

The  second  term  on  the  left  is  negative  since  6  >  0 
z(M  +  M**^)  -  (2xD  +  p(l-5)B,(l-6)Bi)  -  (l-6)(pB,Bi)  >  0 
Hence  the  first  term  of  (3.9)  is  non-positive  so  we  have 

(3.10)  -cx  +  pb  <  0  . 

Now  from  (2.4)  we  have 

(3.11)  (A  -  6B)x  ^  b  for  some  x  >  0 

and  from  (3.6)  p(A  -  6B)  >0  so  multiplying  (3.11)  by 

0  <  p(A  -  6B)x  £  pb  . 

Therefore,  from  (3.10)  we  have  cx  >  0  but  from  (3. 
Now  let  x  be  any  solution  of  (2.2).  Then  x  +  Xx  ,  1  > 
solution  and 

«(x  +  Xx)  *  <i)(x)  +  Xcx  -  XxDx  -  X  xDx 
■  a)(x)  +  Xcx  from  (3.5) 

which  is  unbounded  in  X  ,  contradicting  (2.5)  ■ 


.  Now 
from  (3.5). 
qz  <  0  or 


p  gives 

6) ,  (B  -  A)x  >  0  . 
0  is  also  a 
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4.  AN  EXAMPLE 

Even  when  a  competitive  program  exists,  the  complementary  pivot 
algorithm  may  fall  to  find  it  as  the  following  example  shows: 


The  pair  x  «  3  ,  p  >  (0,1)  is  competitive  because 

(A  -  B)x  -  ('3)  <  (‘3)  ,  p(A  -  6B)  -  (0,1)(^  ■  -  1  -  c  and 

complementary  slackness  holds. 

The  Initial  tableau  of  the  complementary  pivot  algorithm  is 


m. 


m-, 


e  q 


( 


1 

0 

0 


0 

1 

0 


0  0  26  -  1  -1 
0-10  0 


(3) 


1 


1  0 


0 


-1 


The  first  pivot  gives 


1 

1 

2 


and  since  e^  leaves  the  basis,  column  a  must  enter,  but  for  6  ^ 
there  is  no  feasible  pivot  and  the  algorithm  terminates  without  finding  the 
solution.  Of  course,  in  this  case  our  condition  (2.5)  is  not  satisfied. 

For  *5  >  Y  the  algorithm,  as  proved,  will  find  the  solution.  Letting 
e  *  26  -  1  ,  the  tableau  sequence  is 


njIv- 


11 


e 


®3 


1 

e 

0 

0 


1 


e 


-1 


-1 


1 

0 


12  0 


1 

e 

0 

0 


0 

1 

0 


1 

e 

2 

5 


“l 

e 

®3 


-10  1  E 

0  0  0  -£ 

@  1  0  -2e 


-1 

1 

2 


0 

1 

0 


1 

1 

3 


0  110 
0  0  0  -e 

110  0 


0 

1 

0 


-1 


1 


-2 


3 

1 

1 


Conclusion; 

Under  economically  reasonable  conditions  the  complementary  pivot 
algorithm  will  always  find  an  optimal  stationary  program. 


> 
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